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Abstract 

Let X = {X(t) : t S T} be a non-centered, unit-variance, smooth Gaussian random 
field indexed on some parameter space T, and let Au{X,T) = {t G T : X{t) > u} be the 
excursion set of X exceeding level u. Under certain smoothness and regularity conditions, 
it is shown that, as u —>■ oo, the excursion probability P{supjgy A(<) > m} can be approx¬ 
imated by the expected Euler characteristic of Au{X,T), denoted by E{x(A„(A, T))}, 
such that the error is super-exponentially small. This verifies the expected Euler char¬ 
acteristic heuristic for a large class of non-centered smooth Gaussian random fields and 
provides a much more accurate approximation compared with those existing results by 
the double sum method. The explicit formulae for E{x(Au(X, T))} are also derived for 
two cases: (i) T is a rectangle and X — EA is stationary; (ii) T is an iV-dimensional 
sphere and X — EA is isotropic. 

Key Words: Excursion probability, Gaussian random fields, Euler characteristic, rectangle, sphere, 
super-exponentially small. 
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1 Introduction 

Let X = {X{t) : t G T} be a real-valued Gaussian random field living on some parameter 
space T. The excursion probability P{supjg 2 ’-^(b ^ been extensively studied in the 

literature due to its importance in both theory and applications in many areas. We refer to 
the survey Adler (2000) and monographs Piterbarg (1996a), Adler and Taylor (2007) and 
Azais and Wschebor (2009) for the history, recent developments and related applications 
on this subject. To approximate the excursion probability for high exceeding level u, many 
authors have developed various powerful tools, including the double sum method [Piterbarg 
(1996a)], the tube method [Sun (1993)], the expected Euler characteristic approximation 
[Adler (2000), Taylor and Adler (2003), Taylor et al. (2005), Adler and Taylor (2007)] and 
the Rice method [Azais and Delmas (2002), Azais and Wschebor (2008, 2009)]. 
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In particular, the expected Euler characteristic approximation establishes a very general 
and profound result, building an interesting connection between the excursion probability 
and the geometry of the field. It was first rigorously proved by Taylor et al. (2005) [see 
also Theorem 14.3.3 in Adler and Taylor (2007)], showing that for a centered, unit-variance, 
smooth Gaussian random field, under certain conditions on the regularity of X and topology 
of T, 

P |supX(t) > = E{x(A„(A, T))}(1 + o(e“"“ )) as u ^ oo, (1.1) 

where x(A„(A, T)) is the Euler characteristic of the excursion set A„(A, T) = {t E T : 
X[t) > u} and a > 0 is some constant. This verifies the “Expected Euler Characteristic 
Heuristic” for centered, unit-variance, smooth Gaussian random fields. Similar results can 
also be found in Azai's and Wschebor (2009) where the Rice method was applied. It had 
also been further developed by Cheng and Xiao (2014b) that (jl.ip holds for certain Gaussian 
fields with stationary increments which have nonconstant variances. However, to the best 
of our knowledge, among the existing works on deriving the expected Euler characteristic 
approximation the Gaussian field X is always assumed to be centered. In fact, the 

study of excursion probability for non-centered Gaussian fields is also very valuable since 
the varying mean function plays an important role in many models. Especially, when the 
Gaussian field is non-smooth, several results on the excursion probability have been obtained 
via the double sum method [see, for examples, Piterbarg (1996a), Piterbarg and Stamatovich 
(1998), Husler and Piterbarg (1999)]. 

In this paper, we study the excursion probability Pjsupjgj^ X(t) > u} for non-centered, 
unit variance, smooth [see condition (HI) below] Gaussian random fields. As the first contri¬ 
bution, we obtain in Theorem 13.51 that, in general, the expected Euler characteristic approxi¬ 
mation (|2.6I1 holds for such non-centered Gaussian fields when T C is a compact rectangle. 
It shows that, comparing with the double sum method for non-smooth non-centered Gaus¬ 
sian fields [see Piterbarg and Stamatovich (1998) for example], we are able to obtain a much 
more accurate approximation for the excursion probability of smooth non-centered Gaussian 
fields such that the error is super-exponentially small. This is because the expected Euler 
characteristic approximation takes into account the effect of X over the boundary of T, which 
is ignored in the double sum method. By similar arguments in Azais and Delmas (2002), 
such approximation can also be easily extended to the cases when T C is a compact and 
convex set with smooth boundary or a compact and smooth manifold without boundary, see 
Theorem 12.71 

To apply the approximation in practice, one needs to find an explicit formula for the 
expected Euler characteristic E{x(A„(X, T))}. Under the assumption of centered Gaussian 
fields, Taylor and Adler (2003) showed a very nice formula for E{x(A„(X, T))} [see also Adler 
and Taylor (2007)] , involving the Lipschitz-Killing curvatures of the excursion set Au{X,T). 
However, there is lack of research to evaluate E{y;(Au(A,T))} for non-centered Gaussian 
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fields. We provide here explicit formulae of E{x{Au{X,T))} for two cases of non-centered 
Gaussian fields: (i) T is a rectangle and X — EX is stationary; (ii) T is an iV-dimensional 
sphere and X — EX is isotropic; see respetively Theorems 13.51 and 13.111 The results show 
that, the mean function of the field does make the formula of E{x(^m(X, T))} much more 
complicated than that of the centered field. In real applications, one usually needs to use the 
Laplace method to obtain explicit asymptotics for E{x(X„(X, T))}. 

2 Excursion Probability 

2.1 Gaussian Random Fields on Rectangles 

We first consider the Gaussian field X = {X{t) : t G T}, where T C is a compact 
rectangle. Throughout this paper, unless specified otherwise, X is assumed to be unit- 
variance and T denotes a compact rectangle. For a function /(•) G G^(T), we write = 
fi{t) and Denote by V/(t) and the column vector (/i(t),..., fN{t)Y' 

and the N x N matrix respectively. We shall make use of the following 

smoothness condition (HI) and regularity condition (H2) for approximating the excursion 
probability, and also a weaker regularity condition (H2') for evaluating the expected Euler 
characteristic E{y(Xu(X,T))} [note that (H2) implies (H2')]. 

(HI). X(-) G C'^(r) almost surely and its second derivatives satisfy the uniform mean-square 
Holder condition: there exist constants L > 0 and G (0,1] such that 

¥.{Xij{t)-Xi,{s)f <Ld{t,sf^, VLsGT, i,j = l,...,X, (2.1) 

where d{t, s) is the distance of t and s. 

(H2). For every pair (t, s) G with t s, the Gaussian random vector 

(X(t),VX(t),Xi,(t), X{s),VX{s),Xij{s),l< i<j<N) 
is non-degenerate. 

(H2'). For every t ^T, (X(t), VX(t), Xjj(t), 1 < i < j < N) is non-degenerate. 

We may write T = ~cx) < ai < bi < oo. Following the notation on page 

134 in Adler and Taylor (2007), we shall show that T can be decomposed into several faces 
of lower dimensions, based on which the Euler characteristic of the excursion set can be 
formulated. 

A face J of dimension k is defined by fixing a subset cr{J) C {1,..., X} of size k {If k = 0, 
we have cr(J) = 0 by convention) and a subset e(J) = {ej,j ^ o'(J)} C {0,1}'^“^ of size 
N — k, so that 

J = {t = {ti,... Hn) ^ T : Oj < tj < bj if j G cr( J), 

tj = (1 — £j)aj + £jbj if j ^ 
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Denote by dkT the collection of all /c-dimensional faces in T. Then the interior of T is given 

O 

by T = d^T and the boundary of T is given by dT = Ujgaj.T’J. For J G d^T, denote by 
VX\j{t) and V‘^X\j{t) the column vector {Xi^{t), ...,and the kxk matrix 
{Xran{t))ra,n&a(J), respectively. 

If X{-) G C'^{T) and it is a Morse function a.s. [cf. Definition 9.3.1 in Adler and Taylor 
(2007)], then according to Corollary 9.3.5 or pages 211-212 in Adler and Taylor (2007), the 
Euler characteristic of the excursion set Au{X,T) = {t G T : X{t) >u} is given by 

N k 

x{Au{X,T)) = Y, E (2.2) 

k=0J&dkT i=0 

with 

fkiiJ) ■= #{t G J : X{t) > u, VX|j(t) = 0,index(V^X|j(f)) = i, 
e*Xj{t) > 0 for all j ^ it(J)|, 

where = 2ej — 1 and the index of a matrix is defined as the number of its negative 
eigenvalues. 

For t G J G let 

Aj(t) = (Aij(f))ijg^(j) := (Cov(Ai(t),Xj(t)))ijg^(j) = Cov(VA|j(t), VX|j(t)), 

{Ji,--- ,J^_fc} = {l,--- ,X}V(J), (2.3) 

• • • AjN-k) ^ ■ tjE* > 0,j = Ji, - ■ ■ , JN-k}- 

Since A has unit variance, Cov(X(t), V^X|j(t)) = —Cov(VX|j(t), VA|j(t)) = —Aj(t), which 
is negative definite. Dehne the number of extended outward maxima above level u as 

(J) := G J : X{t) > u,VX\j{t) = 0,mdex{V^X\j{t)) = k, 
s*Xj{t) > 0 for all j ^ ci'(J)} 

= {t G J : X(t) > tt, VA|j(t) = 0, index(V^A|j(t)) = k, 
(Xj,(t),...,Aj^_,(t))GE(J)}. 

By similar arguments in Piterbarg (1996b) or Cheng and Xiao (2014b), we have the following 
bounds for the excursion probability: 

E E E{Mf (J)} > p|supA(t) > u\ 

k=0 JedkT ^ ' 

E (E{Mf(J)}--E{Mf(J)(Mf(J)-l)} (2-4) 

fc=o JedkT A 7 

- E IE{Mf (J)Mf (J')}. 

J^J' 

We call a function h{u) super-exponentially small [when compared with P(supjg' 7 i A(t) > 
n)], if there exists a constant a > 0 such that h{u) = o(e“““^““^/^) as n —)■ oo. The 
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sketch for proving the expected Euler characteristic approximation consists of two steps. 
The first step, which is established in Lemma 12.11 below, is to show that the difference 
between the upper bound in (|2.4p and the expected Euler characteristic K{x{Au{X,T))} is 
super-exponentially small. Then we prove that the upper bound in (12.41) makes the major 
contribution since the last two terms in the lower bound in (12.41) are super-exponentially 
small, see Corollary 12.31 and Lemmas 12.41 and 12.51 below. 

Lemma 2.1 Let X = {X(t) : t G T} be a Gaussian random field satisfying (HI) and (H2'). 
Then for each J G d^T with k > 1, there exists some constant a > 0 such that 

E {Mf (J)} = E (2.5) 

Proof To simplify the notation, without loss of generality, we assume cr(J) = {1,... , fe} 
and that all elements in e(J) are 1, which implies E[J) = Let Pj be the collection 

of all k X k matrices with index i. By the Kac-Rice metatheorem [cf. Theorem 11.2.1 or 
Corollary 11.2.2 in Adler and Taylor (2007)], E{M,f(J)} equals 

^E{|detV2X|j(t)|l|v2X|j(t)G(Dfc}l{X(t)>n}l{(Xfc+i(0,-,Xjv(i))GR^^""}l'^^l'^(^^ = 0} 

^ P\/X\j{t)i0)dt 

p poo poo poo 

= (-1)^ dt dx dyk+i ■■■ dvM 
Jj Ju Jo Jo 

E{deW‘^X\j{t)1^^2x^j(t}eVk}\^{t) = x,Xk+i{t) = yk+i, ■ ■ ■ ,X]\[ = yN,XX\j{t) = 0 } 

X Pxit),Xk+iit),-,XN(t){x,yk+i, - ■ ■ ,yN\'^X^j{t) = 0)pvX|j(l(0)- 

(2.6) 

Since Aj(t) is positive definite for every t G J, there exists a kxk positive definite matrix 
Qt such that QtAj{t)Qt = Ik, where Ik is the kxk identity matrix. We write X‘^Xyj{t) = 

and let a\j{t) = CoY{Xi{t),{QtV‘^X\j{t)Qt)ij) for Z = 1, • • • ,N. By 
the well-known conditional formula for Gaussian variables, 

E{{QtX‘^X\j{t)Qt)ij\X{t) = x,VA|j(t) = 0,Xfc+i(t) = yk+i,... ,XN{t) = yN} 

= {QtV^m\jit)Qth + i-6,j,a}^it),...,af^mCoyiX{t),XXit)))-^ (2.7) 

• (x,0, ...,0,yfc+i,--- ,yN)'^- 
Make change of variables V{t) = (V)j(t))i<i,j<A;j where 

Vij{t) = {QtV‘^X\j{t)Qt)ij - {QN‘^rn\j{t)Qt)ij + x6ij, 

i.e., 

QtV^X\j{t)Qt = V{t) + QtV^m\j{t)Qt - xh- (2.8) 
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Denote the density of 


{{Vij{t))i<i<j<k\X{t) = x,VX\j{t) =0,Xk+i{t) = ,XN{t) = yN) 

t’y = (% : 1 < * < j < It follows from (j2.7h and the 

independence of X{t) and XX{t) that independent of x. Let (vij) be the 

abbreviation of matrix Applying ()2.8p yields 


E{det(QtV^A|j(t)Qt)l{v2X|j(t)GX'fc}l-^(^) = VX|j(t) = 0, 

Afc+i (t) = yk+i , • • • , Xn (t) =yN] 

= E{det(QtV^X|j(t)Qt)l{QjV 2 X|^(t)Qt 6 Dj^}|A(t) = x,XX\j{t) = 0, ^ 2 . 9 ) 

Xk+i{t) = yk+i, ■ ■ ■ ,X]si{t) = y^} 

= / det{{vij) + QtX'^m\j{t)Qt-xlk)ht,y^^^,...,yf^{v)dv. 

J {v:{vij)+QtX^rrnj{t)Qt-xIke:'Dk} 

Since is continuous in t and T is compact, there exists some constant c > 0 

such that the following relation holds for all t G T and x large enough: 

x 


(%•) + QN m\j{t)Qt - xlk G Vk, V ||(%)|| < 


( 2 . 10 ) 


Let 


W{t,x,yk+i, - ■ ■ ,2/Af) 
J\v\{ 


det {{vij) + QN‘^m\j{t)Qt - xlk) ht,yk+i,...,yNi'^) dv. 


' { 1 ;: {vij )+Qt J {t)Qt-xlk ^'Dk} 

Then (|2.9p becomes 

/ det {{vij) + Qt^‘^m\j{t)Qt - xlk) ht^yk+i,...,yNiv) dv - W{t,x,yk+i, - ■ ■ ,2/Af)- 

jRk(k + l)/2 

( 2 . 11 ) 


It follows from (I2.10p that 

roo roo 

I{t,x):= / dyk+i--- dyNPx{t),Xk+i{t),-,Xkf{t)ix,yk+i,--- ,yN\'^Xij{t) = 0 ) 
Jo Jo 


X \W{t,x,yk+i, - ■ ■ ,yN)\ 


< 


roo roo 

/ dyk+i--- dyNPx{t),Xk+i{t),-,Xf,{t)ix,yk+i,--- ,yN\'^Xij{t) = 0) 

Jo Jo 


|det {(vij) + QtV^m|j(t)(5t - xlk) \ ht, 


i/fc+i.-" ,yN 


{v)dv 


< Pxit){x\'^X\j{t) = 0) 


IIKi)ll>f 


|det {(vij) + QN^mij(t)Qt - xlk) \ ft{v)dv, 


where ft{v) is the density of ((Vij(t))i<j<j<fc|X(t) = x, VX|j(t) = 0) and the last inequality 
comes from replacing the integral domain by R^“^. Hence there exists some a > 0 
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such that fj I{t,x)dxdt = o{e as u ^ oo. Plugging this, together with (|2.9p 

and (|2.11l) . into (|2.6p . we see that (J)} becomes 


/* ^OO ^OO POO 

(-1)^ / det(Aj(t))dt / dx dyk+i--- / dyN 

Jj Ju Jo Jo 

E{det(QtV^X|j(t)Qt)l{v2X|j(t)e©fe}l"^(*) = J^k+i{t) = yu+i, • • • , Xn = yN,^X\j{t) = 0} 

Px{t),Xk+i(t),-,XM(t){x,yk+i, - ■ ■ ,yAr|VX|j(t) = 0)pvX|j(t)(0)dt. 


= (-1)' 


P poo poo poo 

dt dx dyk+i ■I dyN 
Jj Ju Jo Jo 


E{detV^X|j(t)|X(t) =x,Xk+i{t) = yk+i,--- ,Xn = yN,XX\j{t) = 0} 


^ Px(t),Xu+i(t),- ,Xr,{t){x,yk+u - ■ ■ ,yN\XX\j{t) = 0)pvX|j(t)(0) + o(e 

= e|(-1)‘^(-1)V.(J)| (1+0(6-““’)). 

where the last line is dne to the Kac-Rice metatheorem and the fact that 
k 

^(-l)*|detV2X|j(t)|l|v2X|j(t)GDd = detV2X|j(t), a.s. 
2 = 0 


— OlU —W 


72) 


□ 


Lemma 12.21 below can be derived from Lemma 4 in Piterbarg (1996b). It will be used to 
show in Corollary 12.31 that the factorial moments of M^[J) are usually super-exponentially 
small. 


Lemma 2.2 Let X = {X{t) : t G T} be a Gaussian random field satisfying (HI) and (H2). 
Then for any e > 0, there exists ei > 0 such that for any J G dkT with k > 1 and u large 
enough, 


E{Mf (J)(Mf (J) - 1)} < e-“'/(2/35+£) ^ g-«V(2-£i)^ 

where fij = sup^gj supg^gfc-i Var(X(t)|VX|j(t), V^X|j(f)e). Here and in the sequel, is 
the unit sphere in . 

Corollary 2.3 Let X = {7f(t) : f G T} be a Gaussian random field satisfying (HI) and 
(H2). Then for all J G dkT, E{M®(J)(M,f (J) — 1)} are super-exponentially small. 

Proof If /c = 0, then M^{J) is either 0 or 1 and hence E{M^(J)(M,f (J) — 1)} = 0. If 
k>l, then, thanks to Lemma 12.21 it snffices to show /3j < 1. Clearly, for every e G and 
f G T, Var(X(t)|VX|j(t), V^X|j(t)e) < 1. On the other hand, 

Var(X(t)|VX|j(t),V2X|j(t)e) = 1 ^ Cov{X{t),X^X\j{t)e) = 0. (2.12) 
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Note that the right hand side of (12.121) is equivalent to Aj{t)e = 0. However, by (H2), 
Aj(t) is positive definite, which implies Aj{t)e 7 ^ 0 for all e E Thus for every e E 

and t & T, Var(X(t)|VX|j(t), V^X|j(t)e) < 1. Combining this with the continuity of 
Var(X(t)|VX|j(t), V^X|j(t)e) in (e,t), we conclude /3j < 1. □ 

By similar arguments for showing Lemma 4.5 in Cheng and Xiao (2014b), one can easily 
obtain that the cross terms K{M^{J)M^(J')} in (|2.4I) are super-exponentially small if J 
and J' are not adjacent. In particular, as the main step therein, Eq. (4.13) is essentially not 
affected by the mean function of the field. We thus have the following result. 

Lemma 2.4 Let X = {X{t) : t gT} be a Gaussian random field satisfying (HI) and (H2). 
Let J and J' he two faces ofT such that their distance is positive, i.e., inf^gj^^gj/ ||s —1|| > do 
for some do > 0. Then K{M^{J)M^{J')} is super-exponentially small. 

Next we turn to the alternative case when J and J' are adjacent. In such case, it is 
more technical to prove that E{M,f (J)M®(J')} is super-exponentially small. To shorten the 
arguments for deriving Lemma 12.51 below, we will quote certain similar results in the proof 
of Theorem 4.8 in Cheng and Xiao (2014b) [or Theorem 4 in Aza'is and Delmas (2002)]. 

Lemma 2.5 Let X = {X(t) : t E T} he a Gaussian random field satisfying (HI) and (H2). 
Let J and J' he two faces ofT such that they are adjacent, i.e., inf^gj^^gj/ jjs —tjj = 0. Then 
E{M^(J)M.jf(J')} is super-exponentially small. 


Proof Let L := J f] J' ^ ih. Without loss of generality, we assume 

a{J) = {l,...,l,l + l,...,k}, a{j') = {l,...,l,k + l,...,k + k' -1}, (2.13) 

where 0<l<k<k'<N and k' > 1. Recall that, if A: = 0, then cr{J) = 0. Under 
assumption ()2.13p . we have J E d^T, J' E d^'T and dim(/) = 1. Assume also that all 
elements in e{J) and £{J') are 1, which implies E{J) = and E{.J') = . 

We first consider the case A: > 1. By the Kac-Rice metatheorem, E{M,f(J)M,f (J')} is 
bounded from above by 

/* /* poo poo poo poo poo poo 

dt ds dx dy dzk+i ■ ■ dzk+k'-i / dwi+i ■■■ dwk 
Jj J J' Ju Ju Jo Jo Jo Jo 

E{|detV2X|j(t)||detV2X|j,(s)||X(t) = x,X{s) = y,XX\j{t) = 0, Alfc+i(t) — 

• • • ,Xk+k'-i{t) = Zfc+fc'-z, VX|j/(s) = 0,X;+i(s) = wi+i,... ,Xk{s) = Wk} (2-14) 

X 0) ^k+i, • • •, Zk+k'-h^, wi+i, ■■■,wk) 

A{t, s) dtds, 

JxJ' 

where pt,sA^ ^k+i, • • •, Zk+k'-i,0, wi+i, ...,Wk) is the density of 

{X{t),X{s),VXAt),Xk+i{t),..., Xk+k'-i{t),VX\j> {s),Xi+i (s),..., Xk{s)) 






evaluated at (x, y, 0, Zk+i,, Zk+k'-h 0, wi+i,. ..,Wk). 

Let {ei,e 2 ,... ,eAr} be the standard orthonormal basis of . For t £ J and s E J', let 
St,s = {s — i)^/||'S — t\\ and let ai{t, s) = (cj, A(t)et^s)- Then 

N N 

^{t)et,s = {ei, A{t)et,s)ei = '^ai{t,s)ei. 
i=l i=l 

Since A(t) are uniformly positive definite for all t £T, there exists some ao > 0 such that 


{et,s, A{t)et,s) > ao 


for all t and s. Let 

Di = {{t, s) £ J X J' : ai{t, s) > Pi}, if I + 1 < i < k, 

Di = {{t, s) £ J X J' : apt, s) < —Pi}, if k + 1 < i < k + k' — I, 

Dq = \{t,s) £ J X J' : '^ai{t,s){ei,et,s) > Po k 

i=l ^ 

where Po, Pi,..., Pk+k'-i are positive constants such that Pq + 'Y}!i=iXi^ ^ “o- Similarly to 
the proof of Theorem 4.8 in Cheng and Xiao (2014b), Do U iJ^Pjpp^Di is a covering of J x J'. 
BymM, 

k+k'-l 


E{M,^(J)Mf (/)} < 


A{t, s) dtds + 

^0 i=l+l 


A{t, s) dtds. 


Di 


It can be shown similarly to the proof of Theorem 4.8 in Cheng and Xiao (2014b) that 
f A{t, s) dtds is super-exponentially small. Next we show that f A{t, s) dtds is super- 
exponentially small for i = I + 1,... ,k. 

It follows from (I2.14p that f A(t, s) dtds is bounded above by 


dtds 


Di 


POO POO 

J dx dwiPx(t),\7X^j{t),Xi{s)yx^j,(s)ix,0,Wi,0) 


X E{|detV2X|j(t)||detV2X|j,(s)||X(t) = x,VX\j{t) = 0,Xi{s) = Wi,VX\j,{s) = 0}. 

(2.15) 

Notice that if a subset B <£ Di satisfies inftgenJ, sesnJ' P — t\\ > rjo for some % > 0, 
then similarly to Lemma [2.41 f A{t, s) dtds is super-exponentially small. Therefore, in the 
arguments below, we only treat the case when t and s are close enough or \\t — s|| —)■ 0. 
There exists some positive constant Ci such that 

Px{t),vx^j{t),Xi{s),vXy,{s){x, 0, ruj, 0) 

= Pvx^j,{s),Xi{t),...,Xi_i{t),Xi+i{t),...,Xi,{t){0\X{t) = x,Xi{s) = Wi,Xi{t) = 0) 

X Px{t){x\^i{s) = Wi,Xi{t) = 0)pXi(s)iwi\Xi{t) = 0)pxi(t){^) (2.16) 

< Ci(detCov(X(t), VX|j(t),Xi(s), VX|j,(s)))-^/2 

{x-^2{t,s)P\ f (Wi -^l(t,s)P 


X exp < — 


2al(t,s) 


exp 


2af(t,s) 
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where 


6(t,s)=IE{Xi(s)|Xi(t) = 0} = 
al{t,s) = Yar:{Xi{s)\Xi{t) = 0) 


mis), 

detCov(Xj(s), Xj(i)) 


C 2 {t,s)=E{X{t)\Xi{s)=Wi,Xiit) = 0}, 
crl{t,s)=Xar:{X{t)\Xiis)=Wi,Xi{t)=0). 


In particular, applying Taylor’s formula to Xi{s) [see Eq. (4.23) in Cheng and Xiao (2014b) 
or Piterbarg (1996b)], one has 


^ 2 {t,s) = E{X{t)\(VXi{t),et,s) = Wi/\\s - t\\ + o{l),Xi{t) = 0}, 


= m(t) + (Cov(X(t),(VXi(t),et,,)),0) I Var((VXHO,eM» ^ 


\i{t) 


= m{t) — 


Wi/\\s-t\\ +o{l) - {Vmi{t),et,s) \ 

-mi{t) ) 

aj(t, s)[u;i/||s - tjj - {Xmi{t),et,s) + o(l)] 


(2.17) 


Xar{{XXiit),et,s)) 


and 


aUt,s) = XaT{X{t)\{VXiit),et,s),Xiit)) + o{l) < 1 - <5o (2.18) 


for some do > 0- 

Also, by similar arguments in the proof of Theorem 4.8 in Cheng and Xiao (2014b), there 
exist positive constants C 2 ,Cs,Ni and N 2 such that 

detCov(VA|j(t),Xi(s),VA|j,(s)) > Csljs - (2.19) 

and 

E{|detV2X|j(t)||detV2X|j,(s)||X(t) = x,VX\j{t) = 0,Xi(s) = VX|j,(s) = 0} 

= E{|detV2X|j(t)||detV2A|j,(s)||X(t) =x,VX|j(t) = 0, 

(VXi(t),et,s) = Wi/\\s - t|| + 0(1), VX|j/(s) = 0} 

+ \wi/\\s-t\\f^+ 1). 

Combining (I2.15p . (I2.16p . (I2.17p . (I2.18p . (I2.19P and (|2.20l) . and making change of variable 
w = Wi/\\s — t||, we obtain that there exist positive constants 0^,0^, Cq and such that 
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f fjj A{t, s) dtds is bounded above by 


C 4 I I dtds\\s — t|| 




dx 


POO 

/ dwi{x^^ + \wi/\\s - + 1) 

Jo 




{wi - mj(s))" 


2alit,s) / 2alit,s) 

P P POO POO 

C 4 dtds||s — i||“^ / dx dw{x^^ + \w\^^ + 1 ) 

J J Di J u Jo 

'^ — Tnli\ ^ 0-id^s)[w-(Vmi{t),et,s)+o{X)] \ ^ 

X rnyi)^ Var((VXi(t),et..)) ) 


X exp 


< 6*4 / / dtds\\s — t 


2 cr|(i,s) j-exp I 2 aj{t,s) 

[x-C5 + Pi{Cew-C7)f 


{w - mi{t,s)f 


\-i 


poo 

/ exp 
Ju 


'Di Ju L 2(1 - 5o) 

X + l)exp I - 


dx 


2 al{t,. 


where CTi(t, s) = cri{t, s)/||s —i||, fhi{t, s) = mi(s)/||s —i|| and we have used the fact ai{t, s) > 
/3i > 0 for the last line. This, in turn, ensures that there exists some (5i G (0, (5o) such that 
for sufficiently large u, 


A(t, s) dtds 


D, 


< C 4 exp < — 


2(1 - di) } J Jf). 

? 2 /-f 2„„2 


Is — t|| ^dtds 




{w - mi{t,s)f 


dw 


< C 4 exp < — 


2(1 - (5i) j J Joi 


s — t 


2 o-f(t,s) 

dtds Texp I - |(^ 7 Vp ^ |^|iV 3 + 

Jo I 2(1 — Oi)J 


Since ||s — t||“* is integrable on J x J', we conclude that f fj^, A{t, s) dtds is finite and super- 
exponentially small. 

It is similar to show that f A{t, s) dtds is super-exponentially small for i = k + 
1,... ,k + k' — 1. The case when k = 0 can also be proved similarly. □ 

Now we can derive our main result of this section. 

Theorem 2.6 Let X = {X{t) : t G T} be a Gaussian random field satisfying (HI) and (H2). 
Then there exists some a > 0 such that the expected Euler characteristic approximation o 
holds. 


Proof The result follows immediately from combining (j2.4l) . Lemma l2.ll Corollary 12.31 
Lemma 12.41 and Lemma 12.51 □ 
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2.2 Gaussian Random Fields on Other Sets 

Adler and Taylor (2007) obtained the expected Euler characteristic approximation (ll.ip for 
centered Gaussian fields living on quite general manifolds. Since the method used in this 
paper is different, and it may require more powerful techniques and careful arguments to 
extend the parameter sets to general manifolds, hence we will not attempt to achieve such 
extension here. However, similarly to Aza’is and Delmas (2002), we can easily extend the 
approximation to the cases of smooth and compact manifolds without boundary or convex 
and compact sets with smooth boundary. 

We first introduce some notation. Let (T, g) be an Wdimensional Riemannian manifold, 
where g is the Riemannian metric, and let / be a real-valued smooth function on T. Then 
the gradient of /, denoted by V/, is the unique continuous vector field on T such that 
g{Vf,^) = if for every vector field The Hessian of /, denoted by V^/, is the double 
differential form defined by V^/(^, (") = Hf — where i and C are vector fields and Vg 

is the Levi-Civita connection of {T,g). To make the notation consistent with the Euclidean 
case, we fix an orthonormal frame {Ei}i<i<N, and let 

V/=(/i,...,/iv) = (Fi/,...,EAr/), 

f {Ei, Ej))i<ij<N. 

Note that if t is a critical point, i.e. V/(t) = 0, then f{Ei, Ej){t) = EiEjf{t), which is 
similar to the Euclidean case. As in the Euclidean space, we denote by d the distance function 
induced by Riemannian metric g, which is also called the geodesic distance on {T,g). 

If X{-) G C‘^{T), where T is a smooth and compact manifold without boundary, and it 
is a Morse function a.s., then according to Corollary 9.3.5 in Adler and Taylor (2007), the 
Euler characteristic of the excursion set Au{X,T) = {t £ T : X{t) >u} is given by 

N 

x{Au{X, T)) = (-1)^ j;(-l)V,(r) (2.21) 

i=0 

with 


^.(T) := #{t G T : X{t) > u,VX{t) = 0,mdex{X^X{t)) = i}. 


If T is a convex and compact sets with smooth boundary, then we have 

N ^ N-l 

x{Au{x,t)) = (-i)^^^(-i)v,(r) +(-1)^-1 ^(-i)v,(ar) 
2=0 2 = 0 

with 


^.(T) := G r : X{t) > u,VX{t) = 0,mde^{V^X{t)) = i}, 
li^idT) := G dT : X{t) > n, VX|ar(i) = 0, index(V2X|ar(t)) = i] ■ 
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By similar arguments for Gaussian fields on rectangles in the previous section, together 
with the projection technique in Azais and Delmas (2002) or the arguments by charts in 
Theorem 12.1.1 in Adler and Taylor (2007), we can obtain the following extension, whose 
proof is omitted here. 

Theorem 2.7 Let X = {X(t) : t & T} be a Gaussian random field satisfying (HI) and (H2), 
where T is a smooth and compact manifold without boundary or a convex and compact set with 
smooth boundary. Then there exists some a > 0 such that the expected Euler characteristic 
approximation dEU) holds. 

3 The Expected Euler Characteristic 

We now turn to computing the expected Euler characteristic E{x(A„(X,T))}. To do this, 
we need some preliminary results on calculations of certain Gaussian matrices. 

3.1 Preliminary Computations on Gaussian Matrices 

The following lemma can be obtained by elementary calculations. See also Lemma 11.6.1 in 
Adler and Taylor (2007) for reference. 

Lemma 3.1 (Wick formula). Let (Zi, Z 2 ,..., Zn) be a centered Gaussian random vector. 
Then for any integer k, 

E{ZiZ 2 ■ ■ ■ Z2k+l} = 0 , 

E{ZiZ 2 ■ ■ ■ Z2k} = • • • E{Zi^i^_.^Zi^i^}, 

where the sum is taken over the {2k)\/{k\2^) different ways of grouping Z\, ..., Z 2 k into k 
pairs. 

Let Ajv = (Ajj)i<jj<Ar and Ejv = (Ejj)i<jj<Ar be two N x N symmetric centered Gaus¬ 
sian matrices satisfying the following properties: 

IE{AjjA/j;} — £fi,j,k,l^ dijSkij 
E{EijEki} = F{i,j,k,l), 

where £ and J- are both symmetric function of i, j, k, I, and 5ij is the Kronecker delta 
function. 

The following result is an extension of Lemma 11.6.2 in Adler and Taylor (2007). It will 
be used for computing the expected Euler characteristic of stationary or isotropic Gaussian 
fields. 
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Lemma 3.2 Let = (-Bij)i<ij<Ar be an N x N real symmetric matrix. Then, under (EJP, 

t_-r\kf2U\\ 

E{det(AA, + B^)}= V , 

(3-2) 

E{det(HAr + B^)} = det(i?Ar), 

where Sj{Bi) denotes the sum of the Q principle minors of order j in Bi, and So{Bi) = 1 
by convention. 

Proof We first consider the case when N is even, say N = 21. Then 

E{det(A 2 z + B21)} = ^r/(p)E{(Aij^ + • • • {^21,121 + 

p 

where p = ,i 2 i) is a permutation of (1,2,-•• ,21), P is the set of the {2l)\ such 

permutations, and r]{p) equals +1 or —1 depending on the order of the permutation p. It 
follows from Lemma l3.II that for k < I, E{Aijj • • • A 2 A;-i,j 2 j,_i} = 0 and 

E{Aii, • • • 

^ 2 k,i 2 k} — ,zi,2,Z2) — 5ii^52i2} X • • • 

Q2k 

X {£{2k — I,i2k-i,‘^k,i2k) — ^2fc-i,i2fc-i'^ 2 ^, 12 *,}) 

where Q 2 k is the set of the {2k)\/{k\2^) ways of grouping (ii,Z 2 , • • • ,* 2 ^) into pairs without 
regard to order, keeping them paired with the first index. Hence 

^ 77(p)E{Aii, • • • A2fe,i2j^2fc+i i*2fe+l • B2l,i2l 
p 

= X] * 2 ) “ ^In<^2i2} ■■■{£i2 k - 1, Z2fc-1,2fc, i2 k) - d2k-i,i2k-i^2k,i2u} 

P \Q 2 fe / 

X H2fc+l,j2fe+i • • • ^2Z,i2i 

= X] • • • (52fc-l,i2fe-l<^2fc,i2fc) 1 B2k+l,i2k+l ■ ■ ■ B2l,i2i 

p \Q 2 fe / 

~ A:!2^ det((i)jj j2A:+i<j,j<2zj) 

where the second equality is due to the fact that all products involving at least one £ term will 
cancel out because of their symmetry property, and the last equality comes from changing 
the order of summation and then noting that the delta functions are nonzero only in those 
permutations in P with (ii, Z 2 , • • •, i 2 k-iP 2 k) = {1,2,... ,2k — 1, 2k). Thus 

E{det(A2i + B 21 )} = ^ S2i-2k{B2i)- 

k=0 
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Similarly, we obtain that when = 2/ + 1 


IE{det(A 2 z+i + B 21 + 1 )} = 


k\ 2 ^ 


S2l+l-2k{B2l+l). 


k=0 

The proof for the first line in (I3.2p is completed. The second line in (13.2p follows similarly. □ 

Let B]\f{ii,... ■ ■ ■ ,in) = {Bijif.)i<j^k<n be the n xn principle submatrix of B^ 

extracted from the ii,... ,in rows and ii,... ,in columns in Bn, where 1 < ii <■■■< in 


Proposition 3.3 Let Ajy and be two N x N symmetric centered Gaussian matrices 
satisfying (EU), and let Bn be an N x N real symmetric matrix. Then for x G M, 


N 


E{det(AAr + Bn - xIn)} = ^ 


'ln/2\ 

E 


(_l)7V-n / L^J ^ + 2k)\ 


(N-nV. . 

n=0 ^ ' \ k=0 


kl 2 ^ 


Sn-2k{BN) X 


N-n 


N 


E{det(Hzv + Bn- xIn)] = 


N-n 


n=0 


where Sj{-) is defined in Lemma 


(3.3) 


Proof Applying the Laplace expansion of the determinant yields 


N 


E{det(A^ + Bn- xIn)] = ^(-1)'^-”E{S„(A^ + .Bzv)}x 
By Lemma 13.21 


N-n 


n=0 


Ln/2J 

E{ 5 „(A^ + = E E 

l<i\<---<in<N k=0 


{-lf{2k)\ 

fc!2^ 


^ k\ 2 ^ ( ^ 1 

k=0 \n—2k) 


Ln/2J 


1_y 


{N -n)\ 


fc =0 


{-lf{N-n + 2k)\ 


Sn- 2 k{BN), 


(3.4) 


Bn—2k {Bn {ill ■■■ finfil: ■■■ fin)) 


where the second equality is due to the observation that the sum on all principle submatrices 
of order n in the first line makes every principal minor of order n—2 k appear [^) („_? 2 A;) / {n- 2 k) 
many times. Plugging this into (13.4p yields the first line in (j3.3h . 

By Lemma 13.21 again, 

E{5'jj(.^jv T Bn)} — ^ ^ det(57v fii^ ■■■ finfii^ ■■■ fin)) — Bu^Bn ). 

l<h<---<in<N 


Plugging this into (j3.4p . with A at being replaced by En, yields the second line in (13.3p . □ 
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Remark 3.4 Let B]\f = 0. Then it can be derived from Proposition 13.31 that 

E{det(A7v + Bn — xIn)} = {—^)^Hn{x), 

coinciding with the result in Corollary 11.6.3 in Adler and Taylor (2007), where Hn{x) is the 
Hermite polynomial of order N. Meanwhile, 

E{det(HAr + Bn — xIn)} = (—l)^x^. 

□ 


3.2 Non-centered Stationary Gaussian Fields on Rectangles 


Let X = {X{t) : t € T} be a Gaussian random field such X{t) = Z{t) + m{t), where Z 
is a centered unit-variance stationary Gaussian random field, m is the mean function of X, 
and as usual, T is a compact rectangle. By classical spectral representation for stationary 
Gaussian fields [cf. Chapter 5 in Adler and Taylor (2007)], the held Z has representation 

Z{t) = [ e*<‘’^>lT((iA) 


and covariance 

C{t) = [ e*<‘’^>z/((iA), 

where IT is a complex-valued Gaussian random measure and ly is the spectral measure sat¬ 
isfying z^(M^) = (7(0) = fj^. We introduce the second-order spectral moments 


HJ 



XiXjiy{dX), 


and for any face J G with k > 1, denote Aj = (Aij)jjgo-(j). Then we have 


Cov(VZ|j(t), VZ|j(t)) = -Gov(Z(t), v2Z|j(t)) = Aj 


and that 

£Q{i,j,k,l) :=E{Zij{t)Zkiit)} = [ XiXjXkXiiy{dX) 

JW.N 

is a symmetric function of i, j, k, 1. 

Recall that for a k x k positive dehnite matrix B, the principal square root of B~^, 
which is usually denoted by is the unique k x k positive dehnite matrix Q such 

that QBQ = Ik. Denote by 'I'(x) the tail probability of a standard Gaussian distribution, 
that is 'I'(x) = (27r)“^/^ e~y^^'^dy. Notice that in (13.5p below, for every {t} G OqT, 
VX{t) G E{{t}) specihes the signs of the partial derivatives Xj{t) (j = 1,... ,N) and, for 
J G dkT with A; > 1, the set {Ji,..., Jn-u} and E{J) are dehned in (12.3p . 
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Theorem 3.5 Let X = {X{t) : t £ T} be a Gaussian random field such that X{t) = 
Z{t) + m{t), where Z is a centered unit-variance stationary Gaussian random field and m is 
the mean function of X. If X satisfies conditions (HI) and (H2'), then 


E{xiAuiX,T))} 

= P{VX(i)££({i)))'f(«-m(i)) + E E 

{t}G9oT k=lj£dkT ^ > 

X 


j^dt j dxexp|-i [(x-m(t))2 + (Vm|j(t))'^Aj^Vm|j(t)] 

X P{(Xj,(t), • • • G E{J)\XX^j{t) = 0} 


(3.5) 


where Sj- 2 i{-) is defined in Lemma HOI and is principal square root of A . 


Proof If J = {f} G OqT, then 

E{/io(J)} = P{X(t) > u,e*Xj{t) > 0 for all 1 < j < iV} 

= P{VX(f) G .B({f})}^'(u - m{t)), 

where the last equality is due to the independence of X{t) and XX(t) for each fixed t. 

Let J G dkT with k > 1 and let T>i be the collection of all x A: matrices with index i. 
Applying the Kac-Rice metatheorem, similarly to the proof of Lemma l2.11 we obtain 
k 

e| ^(-i)v,(j)| 

i=0 ' 

„ k 

= / PvX|j(t)(0)dt^(-l)*P{|detV^X|j(t)|l{v2js:|jp)ei>i} 
d i=0 

= j P\7X^j(t){^)dtE{detV‘^X^j{t)l^x(t)>u}'^{{Xj^(t),-,Xjj^_i^{t))eE{j)}\'^^\j{t) = 0 } 

= (2^)(‘+»/=(det(A,;))V2//‘/ [{:c-m(t))^ + (Vm|,(t)fA7>Vm|,(t)]} 

xP{(Xj,(t),-- - ,Xj^_^{t))£E{J)\VX\j{t) = 0}E{detV^X\j{t)\X{t)=x}, 

(3.7) 

where the last equality is due to the fact that VX(t) is independent of both X{t) and V^X(t) 
for each fixed t. 

Now we turn to computing E{detV^X|j(t)|X(t) = x}. To simplify the notation, let 
Q = Then 

E{Z{t){QV^Z\j{t)Q)ij} = -iQAjQfj = -6ij (3.8) 
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and we can write 


E{det{QV^X\j{t)Q)\X{t) = x} = E{det{QV^ Zij{t)Q + QV^mj{t)Q)\X{t) = x} 

= E{det(A(x) + QV‘^mj{t)Q)}, 

where A(x) = (Aij(x))jjgo-(j) is a Gaussian matrix. Applying the well-known conditional 
formula for Gaussian variables and (|3.8D . we obtain 

E{A,,(x)} = E{(QV2Z|j(t)Q),,|A(t) = x} = -x6^j 


and 

E{[Aij(x) - E{Ay (x)}][Afci(x) - E{Afcz(x)}]} 

= E{{QX^Z\j{t)Q),j{QX^Z\j{t)Q)ki} - 6^j6ki = S{i,j,k,l) - 6,,6ki, 

where £1 is a symmetric function of i,j,k,l. Therefore, 


E{det(gv2A|j(t)Q)|A(t) = x} = E{det(A + QV^mj{t)Q - x4)}, 
where A = (Ajj)j jgo-(j) and Ajj are Gaussian variables satisfying 

E{Aij} = 0, E{AijAki} =£{i,j,k,l) - 5ij6ki. 

It follows from Proposition 13.31 that 

E{det(gv2A|j(t)g)|A(t) = x} 




(-i)fe-. {-iy{k-j + 2i) 


-Sj-2iiQX^mj{t)Q) X 


k-j 


j=0 


\ i=0 


Therefore, 


E{detV2A|j(t)|A(t) = x} = det(Aj)E{det(gv2A|j(t)Q)|A(t) = x} 


Plugging this into (j3.7p . together with (j3.6n and (12.2p . yields the desired result. 


□ 


Corollary 3.6 Let the conditions in Theorem |g.5l hold. Assume additionally that to, an 
interior point in T, is the unique maximum point ofm{t) and that V^m(to) is nondegenerate. 
Then as u ^ oo, 

E{x(A4.Y,r))} = ^^^^<P(„-™(i„))(l + o(l)). (3,9) 

\/det(-V^m(to)) 
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Proof By Theorem 13.51 


E{x(A^{X.T))} = 

X J exp | —i [(x — m(t))^ + (Vm(t))^A“^Vm(t)] I (ii(l + o(l)). 
Applying the Laplace method [see, e.g., Wong (2001)], we obtain that as x —>■ oo, 

J exp 1 “^ [(^ “ + {VA~^Vm{t)] | 

(27r)^/2 


dt 




expi--(a:-m(to)) Ml + o(l)) 


Thus as u —>■ oo, 


E{x{AuiX,T))} = j x^/‘^expi-^{x-m{to)f\dx{l + o{l)) 

V27rvdet(—V^m(to)) Ju I ^ J 

- m{to)){l + o(l)). 


'J^\J det(—V2m(to)) 

Y^det(Aj)ri'^/^ 

^det(-V2m(to)) 


□ 


Remark 3.7 The asymptotic approximation in (13.Op is a special case of Theorem 5 in Piter- 
barg and Stamatovich (1998) when the index a therein equals 2, which implies the Gaussian 
held is smooth. However, in our result, a higher-order approximation is also available by 
applying a higher-order Laplace approximation to E{x(Atj(X,T))} [see, e.g., Wong (2001)]. 
Since the calculation is tedious, it is omitted here. □ 


Corollary 3.8 Let the conditions in Theorem \3.5\ hold. If Z is an isotropic Gaussian random 
field with Var(Zi(t)) = 7 ^, then 


E{x(A4X,r))} 


N 


F{VX{t)eE{{t})}d>{u-m{t)) + Y E 


7 


{tj&doT 
X 


k=l JGdkT 


(27r)(^+^)/2 


j^dt j dxexpj-i [(x-m(t)) 2 -h 7 2||Vm|j(t)||2] 


xFiiXjfit),--- ,Xj^_^it))eE{J)} 


X 


Ti/21 




Proof The result follows immediately from Theorem 13.51 the independence of Xi{t) and 
Xj{t) when i j and that, for J E d^T with fc > 1, Aj = 7 ^/fc , which implies Aj = 7 ^4- 
□ 
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3.3 Non-centered Isotropic Gaussian Fields on Spheres 


Let denote the A^-dimensional unit sphere and let X = {X{t) : t E be a Gaussian 
random field such X{t) = Z{t) + m{t), where Z is a centered unit-variance isotropic Gaussian 
random field on and m is the mean function of X. 

The following theorem by Schoenberg (1942) characterizes the covariances of isotropic 
Gaussian fields on [see also Gneiting (2013)]. 

Theorem 3.9 A continuous function C{-, •) : x —>■ M is the covariance of an isotropie 

Gaussian field on if and only if it has the form 

OO 

C{t,s) = '^anPni{t,s)), t,sGS^, (3.10) 

n=0 

where A = {N — l)/2, an > 0, J2'^=o i^) < oo, and Pn is the ultraspherieal polynomials 

defined by the expansion 

OO 

{l-2rx + r^)-^ = Y,r^Pn{x), xE [-1,1]. 

If X is centered, then it only depends on the covariance function which behaves isotropi¬ 
cally over S^. Therefore, as discussed in Cheng and Xiao (2014a) and Cheng and Schwartz- 
man (2015), we do not need to introduce any specific coordinate on sphere. However, if X 
is non-centered, then then mean function m can be very general and hence it is much more 
convenient to use the usual spherical coordinates for arguments, especially for obtaining ex¬ 
act asymptotics. To achieve this, for t = (fi,--- Pn+i) £ we define the corresponding 
spherical coordinate 9 = (0i,..., 0jv) as follows. 

ti = cos 01 , 
t 2 = sin 01 cos 02 , 

^3 = sin 01 sin 02 cos 03 , 


tw = sin 01 sin 02 • • • sin 0Ar-i cos 0jv, 
tN+i = sin 01 sin 02 • • • sin0iv_i sin0jv, 

where 0 E 0 := [OjTt]-^”^ x [0,27r). We also define the Gaussian random field X = {X{9) : 
9 E 0} by X{9) = X(t) and denote by C the covariance function of X accordingly. Similarly, 
let Z{9) = Z{t) and m(0) = m{t). 

Lemma 13.101 below, characterizing the covariance of (X(0), VX(0), V^X(0)), can be ob¬ 
tained easily by elementary calculations. The proof is omitted here. 


20 




Lemma 3.10 Let X = {X{t) : t E §>^} be a non-centered isotropic Gaussian random field 
with covariance (13.101) and satisfying (HI) and (H2'). Then 


dcje,^) _ d^c{e,ip) ^ 

dOi deid(pjd(pk~‘^ 

dOidifj ddidej 


d^Ci9,g^) _ 

d9id9jdipkdipi 


C” {fiijdki + dikSji + 6 ii 6 jk) + C' 5ij5ku 


where 


^^Pni 


n=l 




and C” = ^ ( -^Pn{x)\x=i 

n=2 ^ ^ 


(3.11) 


Now we can formulate the expected Euler characteristic of non-centered Gaussian fields 
on sphere as follows. 


Theorem 3.11 Let X = {^(t) : t E S^} be a Gaussian random field such that X{t) = 
Z{t) + m{t), where Z is a centered unit-variance isotropic Gaussian random field on with 
covariance (j3.10p andm is the mean function ofX. If X satisfies conditions (HI) and (H2'), 
then 


E{x{A^{X,S^))} 

hfiTyl j da;exp|-i [(x-m(0))2 + (C')"iVm(0))|p 


(27r)(^+^)/2 

N 


- (-l)J M)'(iV - j + 2i)\ 

2^ (N - 11! 1 2^ 

\ i=0 




i!2* 


{C')^-^+\C' - iySj_2i {V^m{9)) I 

(3.12) 


where 0 = [0,7r]'^ ^ x [0,27r), 4){9) = ni^i^(siii 0i)^ *, and C' and Sj- 2 i{-) are defined 
respectively in (j3.1ip and Lemma [Ql 


Proof Since is a smooth and compact manifold without boundary, it follows from (j2.21l) . 
Lemma 13.101 and the Kac-Rice metatheorem that 


E{xiAuiX,§^))} 

= (-1)^ / md9 / dxp^^(,)(O)p^(,)(x)E{detV2x(0)|X(0) = x} 

Jo J u 

= (2,r)(^+b?2(C0^/2 X ^ C'-^\\Vm{9))f] | 

X E{detV2l(0)|X(0) =x]. 


(3.13) 


We only need to compute E{detV^X(0)|X(0) = x}. 
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Case 1: C > 1. By Lemma 13.101 similarly to the proof of Theorem 13.51 we get 
E{detV2x(0)|X(0) = x} = E{det[V2z(0) + V^m{e)]\X{e) = x} 

= (C'2 - C')^/^E{det[(C'2 - + (C'2 - C')-^/^V^m{e)]\X{e) = x} 

= (C'2 - C")^/^E{det[A + (C'2 - C')-^/^X^m{e) - C'{C'^ - C')-^/^xIn]}, 


where A = and Ajj are centered Gaussian variables satisfying 

E{A,,AfcJ = (C'2 - C')-^E{Z,,ie)Zkii0)\Xi0) = x} 

= (C'2 - CT^[C”{6,jSki + d.kSji + 5u5jk) + C'S^jSki - C'HiM 

— ^kl 1 

and T is a symmetric function of i,j,k, I. It then follows from Proposition 13.31 that 
E{detV2A(0)|A(0) = x} 

V^m(0) 




i=o 

C'x 


VC'2 - C 


(3.14) 


,VC'2 -cc 

Case 2: < 1. It follows from similar discussions in the previous case and a slightly 

revised version of Proposition 13.31 that 

E{detV2A(0)|X(0) = x} 


_ /W (-^I2\NI2S^ ^ “ J + 2*)! 

C'x 


\ i=0 


'j-2i 


V2m(0) 
VC' - C'2 


(3.15) 


,VC'-C'2. 

Case 3: C = 1. By Lemma 13.101 again. 

fK{detV^X{0)\X{0) = x} = E{det(V2z(0) + V^m{0))\X{0) = x} 

= E{det(“ + V‘^m{0) - xIn)}, 

where E = (Hjj)i<jj< 7 v and Eij are centered Gaussian variables satisfying 


E{E,jEki} = E{Z,ji0)Zkii0)\Xi0) =x} = C''{5iVki + + &iVjk) = HiJVJ), 

and J-” is a symmetric function of It then follows from Proposition 13.31 that 

N 

E{detV2A(0)|A(0) = x} = ^(-l)^-^S’j {V^m{0)) x^-\ (3.16) 

i=o 

Plugging respectively (|3.14p . (13.151) and (13.161) into (|3.13l) . we see that the expected Euler 
characteristic for all three cases above can be formulated by the same expression (j3.12p . □ 
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Remark 3.12 Let m{6) — 0. Let Uj = be the spherical area of the j-dimensional 

unit sphere. Notice that Hermite polynomials have the following properties: 


roo 

/ Hn{x)e~'^^/‘^dx = Hn-i{u)e 

J U 


-V? 12 


Ln/2J 

x” = n\ 

k=0 


kl2’^{n-2k) 


'Hn—2k {x ^, 


where n > 0 and R_i(x) = \/^\l'(x)e*^/^. Applying Theorem 13.111 together with the 
properties above and certain combinatorial tricks, we obtain 

Liv/ 2 J 




-itV2 


N 


j=0 

where poiu) = ^(u), Pjiu) = (27r)“b'i-i)/2^^._^('^jg--!i^/2 £qj. j > £ ^nd 




2 rN\ tjjv if — j is even, 

Vj/ojjv-i ’ 

otherwise 


(for j = 0,1,...,A^) are the Lipschitz-Killing curvatures of [cf. Eq. (6.3.8) in Adler 
and Taylor (2007)]. This coincides with the formula of the expected Euler characteristic 
for centered isotropic Gaussian fields on sphere obtained in Cheng and Xiao (2014a) via a 
geometric approach. However, the result in Cheng and Xiao (2014a) is still more general 
for studying centered isotropic Gaussian fields on sphere since it is also applicable when the 
parameter sets are subsets of S'^. □ 
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